Current-voltage characteristics of semiconductor/ferromagnet junctions in the spin 

blockade regime 
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It was recently predicted [Phys. Rev. B 75, 193301 (2007)] that spin blockade may develop at 
nonmagnetic semiconductor/perfect ferromagnet junctions when the electron flow is directed from 
the semiconductor into the ferromagnet. Here we consider current- voltage characteristics of such 
junctions. By taking into account the contact resistance, we demonstrate a current stabilization 
effect: by increasing the applied voltage the current density through the junction saturates at a 
specific value. The transient behavior of the current density is also investigated. 



There is currently a great deal of interest in 
spin-dependent transport phenomena in semi- 
conductors and their junctions with ferromag- 
neta 1 ! 2 ! 3 ! 4 ! 6 ! 6 !^ 8 ^ 10 ! 11 . 12 ! 13 ! 14 ! 16 . In large part, this 
interest is motivated by the goal to exploit these 
phenomena in new technologies, such as spintronics 
and quantum computation 1 . Recently, some attention 
has been focused on the problem of the extraction of 
spin-polarized electrons from the semiconductor to the 
ferromagne t 11 ' 12 ' 13 ' 14 ! 15 . Despite the apparent similarity 
of spin extraction with spin injection, spin extraction 
shows unique features. In particular, we have recently 
predicted 2 that the spin extraction process at nonmag- 
netic semiconductor/perfect ferromagnet junctions can 
be limited by spin blockade. The physical mechanism for 
spin blockade is the following: the outflow of majority- 
spin electrons from the semiconductor leaves a cloud of 
minority-spin electrons, which limits the majority-spin 
current through the junction. 

In this letter we explore consequences of this phe- 
nomenon that can be easily verified experimentally. In 
particular, we study the current- voltage characteristics of 
such junctions. We do so by considering the conductiv- 
ity of each of its components (semiconductor, ferromag- 
net and their contact). We show that the current flow- 
ing in a circuit involving a semiconductor/ferromagnet 
interface in the spin blockade regime saturates with in- 
creasing applied voltage. Therefore, such an interface 
can be potentially used as a spin-based current stabilizer. 
We also show that in structures with a semiconducting 
region longer that the spin diffusion length, the current 
density saturates to the critical current density j c found 
m Ref. Instead, in junctions with the semiconductor 
region shorter than the spin diffusion length, the asymp- 
totic current value mey be different from j c depending 
on how the semiconductor is connected from the oppo- 
site side of the junction. In particular, if this second 
contact is a good contact with a normal metal, then the 
asymptotic current value is higher than j c . We also con- 
sider transient processes, which, due to the finite response 
time of the spin polarization to the applied voltage, limit 
the speed of operation of such devices. 

The circuit we have in mind is shown schematically in 
the inset of Fig. Q] We consider a voltage source (bat- 
tery) connected to the semiconductor and ferromagnet 



regions of the junction. Assuming that the ferromagnet is 
a good conductor we can neglect the voltage drop across 
it. We also assume a good contact of the voltage source 
with the semiconductor (ohmic or nonlinear contact at 
this junction can be easily incorporated into our model). 
Therefore, there are two components of the total circuit 
where the voltage mainly drops: the semiconductor part, 
and its contact with the ferromagnet. We can then write 
the total applied voltage V as V = V s + V c , where V s 
and V c are voltage drops across the semiconductor region, 
and the contact, respectively. In our model, we consider 
a perfect ferromagnet, such as a half- metal ferromagnet. 
While both spin-up and spin-down electrons arc injected 
from the battery into the semiconductor, only, let say, 
spin-up electrons are extracted from the semiconductor 
into the ferromagnet. 

Spin and charge transport of a non-degenerate electron 
gas in the semiconductor can be conveniently described 
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FIG. 1: (Color online) Current-voltage characteristic of the 
system calculated for several values of the ratio between the 
contact resistance and the semiconductor resistance p^/(p 3 L). 
j c = eNo^/ D/ (2r a f) is the critical current density, — e is the 
electron charge, No is the electron density in the semiconduc- 
tor, D is the diffusion coefficient, and r a / is the spin relaxation 
time. Other symbols are defined in the text. Inset: schematic 
drawing of the circuit. 
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FIG. 2: (Color online) Spin-up density at the junction as a 
function of the applied voltage for several values of p° c /(p s L). 



within the drift-diffusion approximation 2,3 . For simplic- 
ity, we neglect charge accumulation effects as in Ref. 0. 
In the semiconductor region we can then write 



j = <?E = eNo/j, 



Vs 

PsL 1 



(1) 



where j is the current density, a is the conductivity, E 
is the electric field, p is the mobility defined via v 'drift — 
piE, L is the length of the semiconductor, and p s is the 
semiconductor resistivity. Next, we consider the voltage 
drop across the contact. The conductivity of the contact 
is proportional to the density of majority spin electrons 
in the semiconductor near the contact, n^{0). Therefore, 
assuming a linear relationship between the current and 
voltage drop across the contact at a fixed spin-up density, 
n|(0), we write 



2n T (0) y c 
No pi 



(2) 



where p° c is the steady-state contact resistivity at V — > 
(when ra T (0) = N Q /2). Combining Eqs. (HO]) we get 



V = V s + V c 
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N 
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(3) 



Eq. (J3J), which couples V and j, must be supplemented 
by the system of drift-diffusion equations for the semicon- 
ductor region whose solution gives «j(0). This system of 
equations consists of the continuity equations for spin-up 
and spin-down electrons, and the equations for the two 
spin currents: 
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itU) = aE + e-DVn T(; ). 



(5) 



It is assumed that the total electron density in the semi- 
conductor is constant, i.e., n-[(x) + n^(x) = Nq. Cor- 
respondingly, the electric field is homogeneous and cou- 
pled to the total current density as j = epNoEo. The 
boundary conditions are: j-\(0) = j, j|(0) = 0, n^{L) — 
H (L)=N /2._ 

In the following, we will consider separately the two 
cases of long (L ^> l s ) and short (L < l s ) semiconductor 
regions, with l s the spin diffusion length defined below. 

(i) L ^> l s . — In this limit, a steady-state solution of 
Eqs. (|4l5p is known 2 -. The spin densities decay exponen- 
tially from the junction to their bulk values of N /2. The 
decay occurs on the length scale of the up-stream spin- 
diffusion length 2 - 3 l s = 2D/ (jlEq + ^p 2 E^ + AD/T sf ^j . 
The spin-up density at the junction is 2 - 
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Substituting Eq. ([6|) into Eq. ([3|) and introducing the 
dimensionless current density j — j / j c , we get a closed 
equation coupling current density and voltage: 



V 



PsLj c 



= 1 



p s L 1 - 



(7) 



Fig. [T] shows solutions of Eq. ([7]) at different values of 
the ratio of the contact resistance to the resistance of the 
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FIG. 3: (Color online) Current-voltage characteristic of a 
junction with L — 10pm showing the current density satu- 
ration at j > j c . Inset: steady-state current density as a 
function of L at the fixed value of V/(p 3 Lj c ) — 50. These 
plots were obtained using parameter values D — 220cm 2 /s, 



), (4) p — 8500cm 2 /(Vs), No — 5 ■ 10 15 cm~ 3 and r. 



10ns. 



3 



2 4 6 8 10 

Time (ns) 

FIG. 4: (Color online) Transient current (solid line) excited 
by application of step voltages (dashed line). The dotted 
horizontal line corresponding to j = j c is a guide for eye. 
Here we used L — 20/im and p®/(psL) — 1. The rest of 
parameters are as in Fig. [3] 

semiconductor region. All curves saturate at j/j c = 1 
with increasing voltage. The saturation occurs faster in 
systems having smaller contact resistance. In Fig. [2l we 
plot the corresponding spin- up density n^(0). It follows 
from Figs. [Jj and [5] that, for the selected values of parame- 
ters, the current density j is quite close to the critical cur- 
rent density j c at voltages for which 2n|(0)/A< r o ~ 10~ 2 . 
For current stabilization applications, by specifying the 
maximum desired deviation of j from j c , one can obtain 
the minimal voltage V m in required for that deviation us- 
ing Eq. cn^- 

(ii) L < l s . — In this limit, Eqs. (|4l5p . supplemented 
by Eq. ([3]), were solved numerically-^. Starting with 
unpolarized electrons in the semiconductor, we have it- 



erated at each time step Eqs. (|4I5[) with the constrain 
imposed by Eq. ([3]). 

In this regime, the current- voltage characteristics have 
a similar saturation behavior as in the case L 3> l s - How- 
ever, the asymptotic values of the current density (t — > 
oo, V —> oo) are higher than j c (see Fig. [3J. This is due 
to the boundary conditions n-f(L) = ni(L) = No/2. Such 
boundary condition describes a perfect contact of the 
semiconductor with a large reservoir of spin-unpolarized 
electrons. These spin-unpolarized electrons facilitate dif- 
fusion of electrons from the contact region, reducing the 
level of spin polarization near the contact and thus in- 
creasing the current density at which spin blockade oc- 
curs. We plot the current density as a function of L in 
the inset of Fig. [3J For the selected set of parameters, 
the current density starts to deviate noticeably from j c 
in structures with L < 20/im. 

Finally, in view of potential applications, it is impor- 
tant to know the transient behavior of the current den- 
sity. To do this, we consider stepwise voltage changes 
as shown in Fig. 2J This illustrative shape of V was 
selected to show the response to both positive and neg- 
ative voltage increments. The resultant current density 
depicted in Fig. 2] exhibits spikes at each change in V. 
The main change in current density occurs during the 
first several hundreds of picoseconds after the voltage is 
applied. Physically, during this time period the electron 
spin polarization adjusts to a new value of the bias. In 
particular, immediately after an increase of V, n^(0) is 
larger than its steady-state value at the same voltage. 
Therefore, accordingly to Eq. ([3]), a positive spike in j 
appears. Similarly, a stepwise decrease of V results in 
a negative spike. We finally note that current density 
spikes can not be fitted by a single exponent. 

This work is partly supported by the NSF Grant No. 
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